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ABSTRACT 

We have mapped [Nell] (12.8 //m) and [SIV] (10.5 /mi) emission from W51 IRS 2 with TEXES 
on Gemini North, and we compare these data to VLA free-free observations and VLT near-infrared 
images. With 0.5" spatial and 4kms _1 spectral resolution we are able to separate the ionized gas 
into several components: an extended H II region on the front surface of the molecular cloud, several 
embedded compact HII regions, and a streamer of high velocity gas. We interpret the high velocity 
streamer as a precessing or fan-like jet, which has emerged from the molecular cloud into an OB star 
cluster where it is being ionized. 

Subject headings: HII regions — ISM: jets and outflows — ISM: infrared 



1. INTRODUCTION 

Most star formation, acro ss the IMF, takes p lace in 
massive star-forming cores (|Lada fc Ladall2003h . The 
disruption of these cores as high mass stars form termi- 
nates the formation of low and intermediate mass stars 
and is a major contributor to the dynamical energy in the 
dense ISM. While much of this disruption results from 
the ionization of HII regions and their subsequent expan- 
sion, it is becoming increasingly evident that dynamically 
significant amounts of supersonic gas are present during 
the v ery early phases of core dissipation (|De Pree et al.l 
2004|. In this paper we present observations of a jet con- 
taining a significant fraction of the kinetic energy of an 
expanding H II region. 

W51 IRS 2, which includes the compact HII re- 
gion W51d and the maser complex W51 NORTH, lies 
within the HII region complex G49.5-0.4 in W51A, 
a very luminous region of massive young stars and 
mole c ular and ionized gas dBiegind 119751 : iGenzel et al.1 
119821: ICarpenter fc Sanders! Il998f) . The distance to 
W51 NORTH has been determined from H2O maser 
proper motion studies to be 6-8 kpc (jGenzel et alJll98lt 
ISchneps et al.f l981: Im ai et aI1l2002D . which is consistent 
with the far kinematic distance for the molecular cloud 
velocity of 57kms _1 . We will use a distance of 7 kpc, at 
which 1 pc subtends 30". 

Radio free-free and reco mbination line o bservations of 
W51 have been made by iMartinl (|1972l ). iGaume et all 
(1993), and lMehringen(ll994h. The 3 .6 cm free-free image 
of IRS 2 shown in iMehringen (|1994 Figure 19) is from 
observations using all four VLA configurations, which 
preserves both large and small scale structure. We show 
a larger portion of that image in Figure 1 (Plate x) and 
the IRS 2 region in Figure 2a. Various names have been 
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given to the emission peaks in IRS 2. We will refer to 
the free-free peaks labeled a-f in Figure 2a as IRS 2a-2f, 
and the compact infrared continuum source as 2i. Our 
IRS 2a and 2b have been called W51di, our 2d is W51d, 
a nd our 2e is W5 1 d2- 

lOkamoto et al.l (|2001h observed infrared fine-structure 
line emission of [Nell] (12.8 /mi), [Arlll] (9.0/mi), and 
[S IV] (10.5 /mi) from the ionized gas, and concluded that 
the gas is ionized by radiation from stars of spectral type 
08-09. At 7 kpc, the ionized gas in IRS 2 re quires a 
Lyma n continuum luminosity of 6.5 x 10 49 sec -1 (|Martml 
11972ft. and the luminosity of IRS 2 is 2 - 4 x 10 6 L Q 
(jErickson fc Tokunagal[l980l : I Jaffe et al.lll985h . 

Near-infrared observations show a cluster of young 
stars coincident with IRS 2 containing 2-8 times 
as many ionizing stars as the trapezium clus- 
ter (iGoldader fc Wvnn- Williams! [i~99l : iHodapp fc Davis! 
2002). Stellar colors indicate an extinction of Ay~25 
mag. An H-band adaptive optics image with 0.05" res- 
olution, with the 3.6 cm contours overlaid, is shown in 
Figure lb. This image was taken with NAOS-Conica 
on the VLT Yepun as part of program 71.C-0344(A), 
and was obtained from the VLT archive. Most of the 
point sources are stars in the IRS 2 cluster. The diffuse 
emission correlates rather well with the 3.6 cm emission, 
although it appears to be cut off on the northern side, 
with the northeastern half of the IRS 2d missing, proba- 
bly because of extinction. The diffuse H-band emission 
probably is dominated by HI Brackett series lines. 

The W51 NORTH maser group includes H 2 0, OH, 
SiO, and NH3 masers. The majority of t he H2O maser 
peak s are in the 'domi nant center' group (|Schneps et alj 
1981: llmai et al.ll2002l ' m' in Fig. 2a), whic h is centered 
on a more compact cluster of SiO masers (|Eisner et al.l 
|2002[) and on the peak of the thermal NH3 emission 
( Ho et al.lll98l . 



2. OBSERVATIONS 

We observed W51 IRS 2 during the Texas Echelon 
Cross Echelle Spectrograph (TEXES) Demonstration 
Science run on the 8-m Gemini North telescope in July 
2006. TEXES is a high resolution cross-dispersed mid- 
infrared (5-25 /im) spectrograp h capable of 0.4 " and 
Skms" 1 resolution on Gemini (|Lacv et alj [2002) . All 
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data from the Science Demonstration run are available 



from the Gemini archive at http://archive.gemini.edu 



Members of the astronomical community are welcome to 
propose to use TEXES. We expect it to be available in 
alternating semesters on Gemini and the IRTF. 

The observations presented here were made with a 
0.6" slit width and 0.15" sampling along the slit. W51 
IRS 2 was mapped at three spectral settings: 10.5 /xm, 
including the 951.43 cm -1 [SIV] line; 12.8 /im, including 
the 780.42cm -1 [Nell] line; and 13.1/im, including the 
763.23 cm -1 [ArV] line. The maps were made by step- 
ping the 4"-long N-S slit in 0.15" steps across 12" in 
R.A. Overlapping scans were made to cover 9" in DEC 
at 10.5 /im and 7" at 12.8 /im. An additional ten points 
were observed at the beginning of each scan to measure 
the sky background. 

A map of the 12.8 /im continuum, from an 0.6 cm -1 in- 
terval next to the [Ne II] line, is shown in Figure 2b. The 
[SIV] and [Nell] line intensities, integrated over Doppler 
ranges of vlsr = -80 to +80 km s , with continuum sub- 
tracted, are shown in Figures 2c and 2d. The [Ne II] im- 
age is similar to the radio free- free image (Fig. 2a), with 
peaks seen at the positions of IRS 2a-c and possibly 2e-f. 
A bright peak is seen on the south edge of IRS 2d, but 
the northeast side of 2d is weak or absent in all of our 
maps. The [S IV] map is somewhat different, with IRS 2a 
and 2b absent and 2d even more cut off. With these 
bright sources suppressed, 2c and the extended emission 
covering most of the region mapped are more prominent 
in [SIV]. The [ArV] line was undetected, and must be 
fainter than [Nell] by a factor of at least 30. 

Our continuum ma ps are consistent with those of 
lOkamoto et al.l (|2001l ). who derive a silicate feature op- 
tical depth of T9.7 w 2 toward most of the ionized gas, 
implying a reddening between our lines of E10. 5-12.8 ~ !■ 
Given the rather uniform reddening, except toward IRS 
2i, where it is much larger, the differences between the 
[Ne II] and [S IV] maps must be due primarily to exci- 
tation or density differences. Of these two lines, [Nell] 
normally dominates in gas excited by stars of spectral 
type 08 and later, and [SIV] dominates in gas excited 
by earlier type stars. However, [SIV] has a critical den- 
sity for thermalization more than an order of magnitude 
less than [Nell], decreasing the [S IV]/ [Nell] ratio in gas 
with n e > 10 4 cm -3 . IRS 2a and 2b, and the south edge 
of 2d, with [SIV] /[Nell] < 1/3, are of relatively low ex- 
citation (or high density), whereas 2c and the extended 
emission have [SIV] /[Nell] m 3. 

The fits format data cubes are available from the Gem- 
ini archive. We show several position-velocity diagrams 
from the [S IV] cube in Figure 3 and broad channel maps 
in Figure 4. Movies of the P-V diagrams and maps are 
available on-line. The locations of the P-V cut lines are 
shown in Figure 2c. The Doppler shift of the positive 
velocity emission is close to that of the IRS 2 molecular 
cloud. There is a trend that the positive velocity emis- 
sion toward center of the map is more redshifted than 
that near the edges. The emission near the center is red- 
shifted from the ambient molecular gas by a few kms -1 . 
That near the edges is blueshifted by ^20kms -1 . The 
most striking aspect of the P-V diagrams is the wide 
range in Doppler shift seen near and south of IRS 2c. 
The 2c emission ridge is blueshifted from the a mbient gas 
by ~110kms -1 . This emission was seen by iMehringerl 



(|1994h in the H92a line, but was not resolved spatially, 
and the bridge between the blueshifted and ambient 
emission was not apparent. This 'bridge' emission makes 
it clear that the blueshifted gas is not just superimposed 
along the line of sight, but is interacting with the ambi- 
ent gas. This is especially clear in P-V cut c, in which the 
Doppler shift of the ambient gas shifts to the blue where 
it connects into the bridge gas, and there is a gap in the 
ambient gas in the region where the bridge is strongest. 

3. INTERPRETATION 

Gas at Doppler shifts near that of the molecular cloud 
is seen throughout the region observed. Given the pres- 
ence of a cluster of massive young stars centered on this 
region, we identify the exteded gas as the surface of the 
cloud which is ionized by stellar UV. The moderate ex- 
tinction to the gas and stars indicates that they lie near 
the front side of the cloud. The pattern of Doppler shifts 
ranging from just redward of the ambient cloud veloc- 
ity near the center of the region and blueward toward 
the edges is consistent with a flow along a b owl-shaped 
front surface of the cloud. IZhu et ail (|2005l ) show that 
the motion of the ionized gas in several compact H II re- 
gions can be explained in this way, and they suggest that 
such a flow can result, at least in part, from the pressure 
gradient in ionized gas compressed by the ram pressure 
of stellar winds. The [SIV] /[Nell] ratio indicates that 
the gas is ionized by stars of spectral type 08 or earlier. 

The large reddening of IRS2i suggests that it lies 
within the molecular cloud behind the ionized gas. In 
addition, 2a and 2b must be separate from the rest of 
the ionized gas, given their relatively low excitation, in- 
dicating ionization by ^09 stars. However, this gas must 
not be very deeply embedded and could even be in a fore- 
ground cloud. Several of the brightest stars in Figure lb 
lie near ionized gas peaks, suggesting that the star cluster 
is still partially embedded. 

The explanation for the highly blueshifted gas in IRS 2c 
is not so clear, but the presence of 'bridge' gas at in- 
termediate velocities connecting the ambient velocity 
and blueshifted emission in Figure 3 suggests that the 
blueshifted gas is in some sort of a jet that passed through 
the ionized surface of the molecular cloud. The jet must 
be neutral (to be unseen in the ionic lines) while inside 
of the cloud and probably is ionized by stellar UV when 
it emerges from the cloud. The presence of high velocity 
H2O and SiO maser peaks 1" south of 2c (near the Fig- 
ure 3b cut line) suggests that the jet is a part of the same 
outflow as that in which the masers form. If this is the 
case, the jet must have precessed by ~90° since the gas 
we see was emitted, since the masers are moving near the 
plane of the sky in an approximately E-W bipolar pat- 
tern, whereas the ionized jet is to the north of the maser 
center and has a substantial velocity component toward 
us. However, the ionized jet may not be the same as the 
one responsible for the 'dominant center' maser peaks. 
There are several other groups of masers in IRS 2, which 
suggests that there are several outflow centers. We may 
be able to localize better the origin of the outflow we see 
from the information in the [S IV] data cube. Just south 
of 2c there is a gap in the cloud surface gas shown in 
Figure 4a, which could be caused by the passage of the 
jet through the surface H II region. If it is, the source of 
the jet must lie farther to the south, and so to the south 
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of the maser center. Observations of a neutral jet tracer, 
such as shocked H2, could localize the source and test 
this interpretation. 

The E-W elongation of the blueshifted emission could 
be a result of the precession of a narrowly collimated jet, 
or the gas could have been emitted in a fan. In either 
case, the relatively small variation in the Doppler shift at 
the peak of the emission along the blueshifted jet suggests 
that the angle of the jet to our line of sight does not vary 
by a large amount along its length. This indicates that 
the distance along our line of sight to the source of the jet 
is at least several times greater than the 0.13 pc length of 
the blueshifted emission. We can estimate the time since 
the jet was emitted from the ratio of the distance to the 
ionized jet from its source to the observed velocity. This 
gives a time of >1200 yr. 

The extent of the intermediate velocity 'bridge' gas in 
Figures 3 and 4 is something of a puzzle. The narrowness 
of the blueshifted ridge indicates that the jet responsi- 
ble for it is well collimated, at least in declination. But 
the extent of the bridge emission is almost as great in 
declination as it is in right ascension. Two explanations 
seem plausible. As the surface of the jet is photoion- 
ized, it would be expected to expand, and it might slow 
down by interaction with the surrounding medium. Al- 
ternatively, the jet might carry along with it material 
from the HII layer on the surface of the molecular cloud. 
IWilliamsl (|2004l ) has suggested that a turbulent magne- 
tized medium should act as a visco-elastic fluid. It would 
be interesting to ask whether this effect could be strong 
enough to account for the acceleration of the bridge gas. 

We can estimate the mass of the ionized gas in 
the blueshifted emis sion from the ionic line fluxes. 
lOkumura et al.1 (|2001h derived an electron density in the 
ionized gas in IRS 2 of 10 5 — 10 6 cm" 3 from the ratio 
of near- infrared [Felll] lines. At that density, the [SIV] 
line should be thermalized, making its flux proportional 
to mass. We then calculate the mass of ionized gas in 
the jet to be 3 x 10~ 3 M© assuming all of the S is S +3 , 
and S/H is solar. Since some of the gas may be in other 
ionization states (in fact the core of the jet most likely 
is neutral, at least when it emerges from the molecular 
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cloud), this mass estimate is a lower limit. Assuming the 
ionized part of the jet includes gas emitted over 500 yr, 
we obtain a bipolar mass flow of > 1.2 x 10~ 5 M Q yr -1 
and an energy flow of >10L Q . For comparison, molecu- 
lar outflows from m assive young stars rang e upward from 
3 x 10" 5 Mgyr 1 ijOarav fc Lizandll999[ ). 

The high velocity streamer in W51 is similar to the ir - 
radi ated Herbig-Haro j ets se en bv lReipurth et al.1 (|1998l ) 
and ICernicharo et al.l |l998), but has an emission mea- 
sure a factor of 10 5 great er than those of the irra diated 
jets in the Orion nebula (jBallv fc Reipurth1l200lD and a 
mass flow rate a factor of at least 10 4 greater than those 
j ets. It more closely resembles the H H 80-81 bipolar jet 
(|Marti. Rodriguez fc Reipurth! fl993T) . which comes from 
a massive star. The HH 80-81 jet has much greater ex- 
tent than that in W51, but is not as luminous in ionic 
emission, likely due to the lack of as strong of a source 
of ionizing radiation. It is highly collimated and shows 
some evidence of precession, although not as much as we 
suggest for the W51 jet. 
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Fig. 1. — a) (background image) 3.6cm free-free continuum emission from W51 A (IRS 1 and IRS 2). IRS2 is near the top of the 
image. The image was made from the VLA ABCD configuration data set described in Mehringcr (1994). The grayscale uses a square-root 
stretch. The image extent is 112" X 112". b) (inset) H-band continuum emission from stars and ionized gas in IRS 2, with 3.6cm contours 
superimposed. The H-band data were obtained from the ESO VLT archive. The grayscale and contours use a square-root stretch. The 
image extent is 15" X 15". 
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Fig. 2. — Radio free-free, infrared continuum, and infrared line maps of W51 IRS 2. Axes are labeled in arcseconds from the infrared 
continuum source IRS 2i at 19 h 23 m 40.10 s 14°31'06.0" (J2000). a) 3.6 cm emission from Mehringer (1994), with various HII and maser 
peaks marked. Contours are at 1,2, ...8,10,. ..26 X 0.0025 Jybeam" 1 or 7.5 X 10 9 Jysr" 1 . HII peaks marked a-f are referred to in the 
text as IRS 2a-f. The curve through c shows the region of blue-shifted emission. IRS 2i is the infrared continuum source, and m is the 
maser center, b) 12.8/jm continuum emission. Contours are at 1,2,...10,12,...20,24,...44 X 0.014 erg/s cm 2 sr cm - 1 . c) [SIV] (10.5/^m) line 
emission. Contours are at 1,2,. ..15 X 0.009 erg/s cm 2 sr. Lines show locations of P-V cuts in Figure 3. d) [Nell] (12.8 fim) line emission. 
Contours are at 1,2,. ..10,12,. ..20,24,. ..48 X 0.023 erg/s cm 2 sr. 
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Fig. 3. — Position- velocity cuts through the [SIV] data cube. Horizontal axes are labeled in arcseconds north from IRS 2i. The contours 
are spaced quadratically. The ambient molecular cloud is near vlsr = 57kms _1 . Locations of cut lines are shown in Figure 2c. The online 
movie shows the data cube as a sequence of P-V diagrams. The movie frames are labeled by the R.A. offset from IRS 2i. 
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Fig. 4. — Images of [S IV] emission in three velocity intervals. The brightness stretch is linear and the same for all images. Annotations 
are as in Figure 2. a) Ionized gas near the molecular cloud velocity. We interpret the dark region at bottom center as the region where 
the cloud surface gas has been disrupted by the jet. b) The 'bridge' gas that we interpret as having been accelerated by the jet. c) 
Blue-shifted emission from the jet which has been ionized by stellar UV radiation. The online movie show the data cube as a sequence of 
180 narrow-band images, running from vlsr = +85kms _1 to -82kms~ 1 



